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ABSTRACT 

We present new models for low-mass stars down to the hydrogen-burning limit that consistently couple atmosphere and 
interior structures, thereby superseding the widely used BCAH98 models. The new models include updated molecular 
linelists and solar abundances, as well as atmospheric convection parameters calibrated on 2D/3D radiative hydrody¬ 
namics simulations. Comparison of these models with observations in various colour-magnitude diagrams for various 
ages shows significant improvement over previous generations of models. The new models can solve flaws that are 
present in the previous ones, such as the prediction of optical colours that are too blue compared to M dwarf obser¬ 
vations. They can also reproduce the four components of the young quadruple system LkCa 3 in a colour-magnitude 
diagram with one single isochrone, in contrast to any presently existing model. In this paper we also highlight the need 
for consistency when comparing models and observations, with the necessity of using evolutionary models and colours 
based on the same atmospheric structures. 
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1. Introduction 


In 

els 


1998, our team released a set 
for low-mass stars (Baraffe et al. 


of evolution ary mod- 
hereafter 


1998 


BCAH98) based on the so-called JNextGen atmosphere 


models ( jHauschildt et al. 1999) that marked a new era 
of models that consistently coupled interior and atmo¬ 
sphere structures. These models became very popular be¬ 
cause they successfully reproduce various observational 
constraints, such as mass-luminosity and mass-radius re¬ 
lationships and colour-magnitude diagrams. The models, 
however, had some important shortcomings, such as pre¬ 
dicting optical ( V—I ) colours that are too blue for a given 
magnitude (see j |3.1[ ). Later generations of models included 
improved molecular linelists for various atmospheric ab¬ 
sorbers, such a s the AMES line l ists u s ed in the Dus t y and 
Cond models dChabrier et al. 


2000 Allard et al., 2001 
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(1998). The major changes concern the atmosphere mod- 
els, which provide the outer boundary conditions for the 
interior structure calculation, and the colours and magni¬ 
tudes for a given star mass at any given age. Substantial 
changes have been made since the NextGen atmosphere 
models used in the BCAH98 evolutionary models. A prelim¬ 
inary se t of atmosphere models, referred to as the BT-S ettl 
models (Allard et al. 2012a|b Rajpurohit et al. 2013), in¬ 
clude some of these changes, which are briefly summarised 
below. More recent modifications concerning the treatment 


of convection are described in 12.3 


Bara ffe et al. 2003p , but they still show shortcomings (see 
p.2p . After a long effort to solve these flaws, efforts have 
paid off with the release of current models that supersede 
the BCAH98 models. In this paper, we describe the main 
physical ingredients of the models and compare them to a 
selection of observations that highlight the improvement of 
these new models over previous ones. 


2. Model description 

Evolutionary calculations are based on the same input 
physics descri bing stellar and substellar inter ior structures 
as are used in Chabrier & Baraffe (19971 and Baraffe et al. 


2.1. Molecular linelists and cloud formation 

Line opacities for several important molecules have been 
up dat ed, notably the water linelist from Barber et al. 
(200 6|, m etal hydrides such as CaH, FeH, CrH, TiH from 
Bernath (20061, vanadium oxide from Plez (200 4 , priv . 
comm.), and carbon dioxide from Tashkun et al. (2004). 
For TiO, t he present s et of atmosphere models uses the 
linelist from Plez (1998). This list is not as c omple te at high 


energies as the AMES 

linelist 

(Schwenke 

1998) adopted 

in 

Allard et al. 

(2001 

2012a 

, with on 

y 11x10° lines 

compared to the 

l&Tx 

10° of 

Schwenke 

(1998 

). But the 


and thus generally improves the optical colours (Fig. [4]). 
Obviously, the field is still in need of a new, complete, 
and accurate theoretical TiO linelist to allow quantita¬ 
tive high-resolution spectroscopic analysis of this impor¬ 
tant molecule, as recently pointed by the high-resolution 
transmission spectrum analysis of a transiting exoplanet 


(Hoeijmakers et al. 20141. 
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Condensation of up to 200 types of liquids and solids 
is included in the atmospheric equation of state. The for¬ 
mation and sedimentation of clouds and depletion of con¬ 
densible mat erial is treate d in the self-consistent timescale 
approach of Allard et al. (2012a). They used a monodis- 
perse grain size distribution for each cloud layer, but cur¬ 
rent version of the models now considers a log-normal distri¬ 
bution over 1.8 decades in size in 12 bins with the altitude- 
depend ent size determ ined by the same timescale formalism 
as in Allard et al. (2012a). The set of 60 grain types included 
in calculating the cloud opacity has also been updated with 
the optical data for several low-temperature condensates, 
which become important in the atmospheres of T dwarfs, 
which are not discussed in this paper (Homeier et al., in 
prep). 


2.2. Solar abundances 


Solar abundances have been revised over the past decade 
based on radiation hydrodynamical simulations of the so¬ 
lar photosphere combined with 3D NLTE radiative trans¬ 


fer (Asplund et al. 


2009 Caffau et al. 


2011). This leads 


to a substantial reduction of carbon, nitrogen, and oxy¬ 
gen a bundances compared to the abundances of Grevesse| 


et al. (1993), previously used in the NextGen and AMES- 
Dusty/Gond models. The present m odels adopt the solar 
composition of Asplund et al. (2009) with revisions of the 
elemental abundances of G, N, O, JNe, P, S, K, Fe, Eu, Hf, 


Os, and Th obtained by the CIFIST project (Ca ffau et al. 
2011). This yields, in particular, a slight upward revision of 


carbon, nitrogen, and oxygen abundances and an increased 
total heavy element fraction by mass, with Z = 0.0153 


compared to 0.0122 in Asplund et al. (2005) and 0.0134 in 
Asplund et al. J 2009 Although the individual differences 
with the Asplund et al. (2009) ele mental abundan c es are 
modest, Antia & Basu (2011) and Basu &; Antia (2013) 
have found that the higher overall metal content of the 
CIFIST abundances is in better agreement with helioseis¬ 
mology results for the solar interior. The helium abundance 
adopted in the atmosphere models is fixed to Fj ni = 0.271, 
which is representative of the initial sola r heliu m abundance 
(Serenelli & Basu 2010 Basu & Antia 2013). 


2.3. Treatment of convection 

Radiative-convective equilibrium in the PHOENIX atmo¬ 
sphere model s is calcula ted usin g the mixin g length theory 
(or MLT, see Kippenhahn & Weigert 1990). The main free 
parameter within this framework is the mixing length Z m i x , 
expressed in terms of the pressure scale height Up and Z m ; x 
describes the efficiency of the convective energy transport in 
terms of the superadiabaticity of the temperature gradient. 
All previous PHOENIX models have used a constant value for 
Z mix = ZjUjg): Z(A“ = Up in t he NextGen and Con d/Dusty 
models and ZjA™ = 2 x Up in Allard et al. (2012a). 


In the present work, the determination of Z'r'-™ is inspired 


by the work of Ludwig et al. (1999 2002). It is based on 
comparisons between ID MLT models and 2D/3D RHD 
simulations that cover main sequence and pre-MS models 


down to th e hydrogen-burning limit and below (Freytag 
et al. 20101 2012 1 . Since the RHD models include a basic 
treatment of dust formation, they also allow the calibra¬ 
tion of late M dwarfs, where cloud opacity is becoming rel- 



log P 



log P 


Fig. 1. Comparisons of pressure-temperature atmospheric 
profiles between present ID (solid blue) and CO s BOLD 
RHD (dash-dot black) models, for different effective tem¬ 
peratures and surface gravities, as indicated in the panels 
for each curves. Pressure unit is dyne/cm -2 . Calibration of 
Z m ix based on such comparisons yields Zf r |™ = 1.64 x Up 
for T eS = 2880, log g = 4.0; C" = 1-58 x U P for T eS = 
3200, log 5 = 4.0; l%£ = 1.69 x U P for T eS = 3470, log g = 
5.0; 1 %£ = 1.58 x U P for T eS = 3960, log g = 5.0. The dot 
on each solid curve indicates the onset of convection and 
the triangle indicates the location of optical depth r =100 


(Chabrier & Baraffe 1997) where the transition to the in- 


terior model is performed. 


evant. This calibration yields a value of ZjA™ ~ 1.6 x Up 
for the Sun. The value of ZA™ increases for later type stars, 
up to values £ 2 xUp for the coolest and densest mod¬ 
els (T e s < 3000 K and logg > 4.5). Comparisons between 
ID and RHD pressure-temperature profiles are shown in 
Fig.[l| The adopted value of Z A™, for a given T e g and log g, 
provides the best overall agreement that can be obtained; 
however, the agreement is not perfect all the way from the 
top to the bottom of the atmosphere, reflecting limitations 
in the MLT formalism for properly handling atmospheric 
convection (Homeier et al., in prep). 

2.4. Evolutionary models 

Current models cover the evolution of pre-MS and MS stars 
from 0.07 M 0 to 1.4 M 0 . Calibration of a 1 M 0 star se- 
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T e „ (K) 


Fig. 2. Comparison of present models with the BCAH98 
models for various isochrones, as indicated in the figure 
and inset. Solid (red): present models; Long dash (black): 
BCAH98 models with = Hp; Dotted (black): BCAH98 
with = 1.9 x Hp. 


quence to fit the fundamental parameters of the Sun (Rq, 
Lq ) at its age requires a value for the interior structure of 
the mixing length Z m ; x = Z“( x = 1.6 x Hp and a helium 
abundance Y = 0.28. We used this value of helium for the 
complete grid of evolutionary models. The value adopted in 
the atmosphere models is slightly different (Y = 0.271, see 
§2.2). To appreciate the effect of this inconsistency, we have 
recalculated a full grid of evolutionary models assuming the 
same He abundance Y =0.271 in the interior structure as in 
the atmosphere models. The effect of such a variation in Y 
on T e g for the whole range of masses considered is less than 
1%. The maximum effect on the luminosity is found for the 
highest masses (M Xl 1 Mq) and reaches at most 10%. For 
a 1 Mq star, the effect on the luminosity is less than 5% 
for ages ,$ 2 Gyr, and up to 10% for ages £ 2 Gyr. We 
have also computed a few atmosphere models again with 
Y = 0.28 for a range of T e g between 2600 K and 6000 K 
and of gravities logg between 3.5 and 5. The effect of such 
a variation in Y is negligible on the atmospheric thermal 
profiles (less than a 1% effect on the pressure at a given 
temperature) and on the photometry (differences less than 
0.01 mag for all colours explored). 

We fixed /)"f x = 1.6 x Hp for the complete grid of evolu¬ 
tionary models for the interior structure. A fully consistent 
approach would be to adopt the value obtained for the cor¬ 
responding atmosphere model, i.e Z™; x = /A™. We checked 
that this is unnecessary. Indeed, values of /A™ only start to 
depart significantly from 1.6 x Hp for the coolest and dens¬ 
est models ( T e g < 3000 K and logg > 4.5), reaching val¬ 


ues up to 2xHp or more. This concerns masses < 0.2 Mq, 
where a variation of 1™\ X between 1.6xHp and 2xHp has 
an insignificant effect on the evolution. 

As illustrated in Fig.[2j the new models predict different 
positions of isochrones than do the BCAH98 models. Below 
T e g ,$ 2800 K, dust starts to form in the upper atmosphere 
and affects atmospheric thermal profiles and spectra. In the 
current grid of evolutionary models with T e g > 2000 K, the 
effects of cloud opacity only gradually become visible for 
the latest M and early L dwarfs (T e g 4 $ 2500 K). They are 
becoming crucial for models of L and T dwarfs, which will 
be examined in a forthcoming study. 

Theoretical mass-radius relationships for 0.1 Gyr and 
1 Gyr ar e shown in Fig. [3| for differ ent sets of models. As 
shown in Chabrier & Baraffe (19971, the mass-radius rela¬ 
tionship for low-mass stars is essentially fixed by the equa¬ 
tion of state, with a very small dependence on the atmo¬ 
sphere treatment. The differences between present models 
and the BCAH98 models essentially stem from the different 
values of the interior mixing length l^ x used. As shown in 
Fig-! the effect of the interior He abundance Yi nt slightly 
affects the mass-radius relationship for masses ^ 1 Mq. 



0.6 0.8 1 1.2 1.4 

M/M 0 


Fig. 3. Comparison of theoretical mass-radius relationships 
for ages of 0.1 Gyr (upper panel) and 1 Gyr (lower panel) 
and based on different sets of models: Present models with 
Fi n t = 0.28 (solid red) and YS nt = 0.271 (short-dash red); 
BCAH98 models with Z(A X = Hp (long-dash black) and 
^mix = 1.9 x Hp (dot black). 


3. Comparison with observations 

In this section, we highlight the significant improvement 
provided by this new set of models, which solve some of 
the problems of the BCAH98 models, through a selection 
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(V-I) 


Fig. 4. Comparison of 1 Gyr isochrones with disk dwarfs 
in an optical CMD. Present models based on the Plez TiO 
(solid red) and the AMES TiO (dash red) linelists are 
shown, along with the BCAH98 models (long dash black). 
The magenta open circles are d isk objects with parallax 
(Monet et al. 1992 Leggett 1992) and the blue squares are 
fr om the recent 5 pc and 10 pc solar neighbourhood samples 
of Cantrell et al. (2013) with more accurate parallaxes and 
cleaned of binary systems. The numbers close to the red 
open circles on the solid curve give T e g- and mass (in M 0 
in the brackets) for selected models. 


of observational tests. Colour-magnitude diagrams (CMDs) 
provide the best validity test, as they test the consistency 
between interior and atmosphere structures. 

3.1. Optical colours 

As mentioned in the introduction, a well-known flaw of the 
BCAH98 models is the prediction of colours that are too 
blue (V—I). Figure El compares models at 1 Gyr with solar- 
neighbourhood and disk objects with parallax (Monet et al. 


1992] Leggett 1992 


2013). The comparison 


Cantrell et al. 

shows that this persistent problem of the models seems to 
be solved by several improvements. The most important one 


is the use of the TiO linelist from Plez (1998). The impact 
of the TiO linelist has already been noticed with the gen¬ 
eration of models following the NextGen ones, such as the 


Dusty models (Chabrier et al. 2000 Allard et al., 2001 


The former ones use the linelist by Jorgensen ( 1994|), whi 


the latter ones include the AMES TiO linelist of ISchwenkel 


(1998). The improvement in (V— I) colours of models that 
indue e the AMES_TiO linelist of Schwenke (1998) was sig¬ 
nificant (see Fig. El), but still unable to yield satisfactory 



Fig. 5. Comparison of 1 Gyr isochrones with disk dwarfs in 
near-IR CMDs in CIT filters. Solid (red): new models; Long 
dash (black): BCAH98 models with Z™j x = Hp. Symbols are 
the same as in Fig. [4] 


agreement with observations. Indeed, an offset in (V — I) 
colours of ~ 0.4 mag remained between T e g ~ 3600 K 

Chabrier] 


see 


(My ~ 10) and T eff ~ 2300 K (My ~ 19,_ 

et al. 2000|). Im provement due to the implementation of 
the |Plez| (1998) TiO linelist in the new models is illus¬ 
trated in Fig7]4 by the comparison between present models 
(TiOPl) and the same models based on the AMES TiO 
linelist (TiOAM). 

Another positive effect stems from the change in solar 
abundances, with a decrease in oxygen abundance by ~22% 


between the 

Grevesse et al. 

(1993[), used in previous sets of 

models, and the presently used 

Caffau et al. ( 

2011 ) abun- 

dances. As a consequence of re 

atively less O, the flux in 


the IR generally increases owing to weaker water absorp¬ 
tion. Because of flux redistribution, the flux in the F-band 
decreases relative to the flux at longer wavelengths. This 
is an important effect, illustrating how water abundance 
affects optical colours due to flux redistribution. 

3.2. Near-infrared colours 

The new models yield also a significantly better match 
with observations for relatively old objects (age ^ 1 Gyr, 
log g 4.5) in near-IR CMDs, as illustrated in Fig. [H] This 
stems fromthe use of th e more complete and accurate water 
linelist of Barber et al. (2006), compared to previously used 
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ones, 


nam ely Miller et al. (1994) in NextGen and Partridge 
& Schwenke ( 1997| ) in Dusty/Cond. 

To illustrate the impact of the new models at low 


gravities and in a similar range of effective temperatures 
(T e ff > 2000K), we also compared them to data in the 
young cluster a Orionis (Pena Ramirez et al. 2012), with 


an age of a few Myr (log g ^ 4). The results are shown in 
Fig- ® in various CMDs. Fluxes from the models have been 
transformed in the VISTA filters to be consistent with the 


data o f|Pena Ramirez et al. (2012). Comparison of models 
in CMDs" (or luminosity-2 eg aiagrams) with observations 
in young clusters and star formation regions should always 
be taken with caution owing to the observed luminosity 
sprea d and the effect of in i tial conditio ns and accretion his¬ 
tory (iBaraffe et al. 2002112009112012 1. The comparison in 


Fig-® however, is illustrative of the overall improvement in 
ZJHK filters of the new models compared to the BCAH98 
ones. The latter indeed provide poor agreement with ob¬ 
servations at young ages in near-IR CMDs, as illustrated 
in Fig. ® For this reason, models had to be compared with 
data in L — T e g or HR diagrams, which implied the use 
of very uncertain bolometric corrections and spectral type- 
re g or colour-T e g scales. The agreement of the new mod¬ 
els with observations, however, deteriorates at the faintest 
magnitudes by up to ~ 0.1 mag in ( J—H ) and ~ 0.4 mag 
in (J — K ). Given photometric uncertainties and possible 
iF-band excess at this age owing to the presence of a cir- 
cumstellar disk, it is premature to draw any conclusion on 
the reason for these discrepancies. 


3.3. Multiple systems 

Stringent tests of evolutionary models are provided by co¬ 
eval multiple systems. The recent discovery of the quadru¬ 
ple pre- main s eq uence syste m LkCa 3 is interesting in this 
respect (Torres et al., 2013). If one admits that the four 
components of LkCa 3 should be coeval and must lie on the 
same isochrone, they make an excellent test of evolutionary 
models. This supposes that early history of accretion has 


not altered their structure (Baraffe et al. 2009 2012) and 


that their position in a CMD or L — T e g diagram can be 
reproduced by non-accreting evolutionary models. 


With this assumption, Torres et al. (2013) find that no 
currently available model is able to reproduce the posi¬ 
tion of the four components with one single isochrone in 
a (V—H) — My CMD. In contrast, our new models do ful¬ 
fil this stringent constraint within the error bars (see Fig. 
[7]). The new models agree with observations for an age of 
1.6 Myr. In comparison, the best-fit 6.3 Myr isochrone of 
the BCAH98 models with Z™; x = Hp is unable to match 
the coolest component LkCa 3 Ab, which can be fitted by 
an isochrone of 2.5 Myr. A similar problem arises with the 
BCAH98 models using = 1.9 x Hp. The same failure 
was reported by Torres et al. ( 2013|) (see th eir Fig. 6 ) with 
the Dartmouth models ( )Dotter et al. 2008). 

As a final note, even if accretion history can have an 
im pact o n the structure of young objects, as suggested by 


Baraffe et al. (2012), this does not necessarily affect all ob- 
jects in the same way. Even for binary or multiple systems, 
each component has its own accretion history, and its struc¬ 
ture may be differently altered by accretion, since the effect 
depends on the total amount of accreted mass (see Baraffe 
et al. 2012 for details). Disentangling the remaining uncer- 



1 2 30.4 0.5 0.6 0.T).6 0.8 1 1.2 

(Z-J) (J-H) (J-K) 


Fig. 6. Comparison of models with observations in a 
Orionis in various CMDs (VISTA filters). Isochrones of 1 
Myr and 10 Myr are displayed for various sets of models: 
Present models (red): solid (1 Myr) and dash-dot (10 Myr). 
BCAH98 models with Q\ x = H P (black): long dashed (1 
Myr) and dot (10 Myr). The data (blue dots) are from 


Pena Rann'rez et al. (2012) using a distance modulus of 
7.73 (i.e. a distance of 352 pc). The hydrogen-burning limit 
(M=0.07 M 0 ) is located at Mj ~ 6 (T e ff ~ 2900K, log <7 ~ 
3.5) at 1 Myr and Mj ~ 8 (T e g ~ 2950K, log g ~ 4.2) at 
10 Myr. 


tainties of models and possible effects of accretion when 
comparing models to young multiple systems is thus not 
straightforward. We hope that many more such multiple 
systems, spanning a wide range of masses and ages, will be 
discovered in the near future to provide more clues to the 
remaining model uncertainties (e.g. convection treatment, 
opacity uncertainty) and/or to the impact of accretion for 
the youngest systems. 

3.4. Tests of stellar evolutionary models 

For the sake of testing evolutionary models against obser¬ 
vations, it is crucial to use the synthetic colours and mag¬ 
nitudes predicted by the same atmosphere models as those 
that provide the outer boundary conditions for the interior 
structure. Using “hybrid” models, based on atmospheric 
thermal profiles for the interior structure from one set of 
atmosphere models and colours/magnitudes from another 
set, yields incorrect results. Figure ® illustrates this point, 
showing in two CMDs for various isochrones a comparison 
between the present set of models, the BCAH98 models, 
and a hybrid model based on BCAH98 evolutionary mod¬ 
els (relying on the NextGen atmosphere models) and the 
colours and magnitudes from current atmosphere models. 
The positions and slopes of the hybrid model isochrones are 
significantly different from both the new and the BCAH98 
models. Using these and similar hybrid models to derive 
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Fig. 7. Comparison of models with the four component s 
of the quadruple system LkCa 3 (Torres et al. 20131. 
The best-fit isochrone is shown for various sets of models. 
Present models for an age of 1.6 Myr (solid red). BCAH98 
models with = Hp for an age of 6.3 Myr (long-dash 
black). BCAH98 models with = 1.9 x Hp for an age of 
4 Myr (dot black). The blue dots are for LkCa 3 A and the 
cyan squares for LkCa 3 B. We adopt a dista nce d = 127 p c 
and extinction E(B — V) = 0.10 as in 
The numbers close to the red open circ. 


Torres et al. (20131. 
es on the soliclcurve 
give T e s and mass (in M 0 in the brackets) for a few models. 



2 3 4 5 6 


(V-H) 



Fig. 8. Comparison of 10 Myr and 1 Gyr isochrones for 
various sets of models in two CMDs. Solid (red): Present 
models. Long dash (black): BCAH98 models with Z("j x = 
Hp. Dot (blue): hybrid models using the BCAH98 structure 
models and the colours/magnitudes of present atmosphere 
models. 


cluster ages and masses or to test them against multiple 
systems yields incorrect results and incorrect conclusions 
on the validity or invalidity of the models. Such compar¬ 
isons, often seen in the literature, are totally meaningless. 


4. Conclusions 

The new set of evolutionary models presented in this pa¬ 
per, consistently coupling internal and atmospheric struc¬ 
tures, show significant improvement over previous gener¬ 
ations. They solve some of the persistent flaws present 
in them, such as predicting optical colours that are too 
blue. Significant improvement in the atmosphere models, 
in terms of updated molecular linelists and revised solar 
abundances, provide a much better match to observations. 
More systematic tests of current models against observa¬ 
tions will help to identify the remaining uncertainties. The 
calibration of the mixing-length parameter, as presented 
in this work, provides overall good agreement between ID 
and multi-D RHD thermal profiles. This exercise highlights 
the poor approximation of using a constant value for (A™ 
for pre-main sequence and main-sequence low-mass stars. 
The agreement, however, is not perfect, depending on the 
effective temperature and the surface gravity. This very 
likely reflects the limitation of the MLT formalism to cor¬ 
rectly capture the thermal properties of convectively unsta¬ 


ble atmospheres over a wide range of parameters. We note 
as well that remaining uncertainties in the multi-D RHD 
simulations are not excluded. Efforts towards developing 
such simulations in cool atmospheres need to be pursued. 
Since models at young ages are particularly sensitive to th e 
treatment of atmospheric convection (Baraffe et al. 2002), 
comparisons between models and observations of objects in 
young clusters (age ^ 10 Myr) will help to improve the 
treatment of convection. 

Following the release of these new models down to the 
hydrogen-burning limit Q we are currently working on their 
extension to L, T, and Y dwarfs, including dust formation 
and settling. This will provide a complete and coherent set 
of models for low-mass stars and brown dwarfs down to 
Jupiter-mass objects. 
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